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The  study  analyzes  the  internal  states  of  a  3-component  system  with  one  active 
element  and  two  spares  in  cold  standby  (pure  replacement  policy  without  repair). 
Ulements  of  the  system  are  assumed  to  have  exponentially  distributed  lifetimes, 
however,  special  attention  is  paid  to  systems  composed  of  components  with  different 
failure  rates.  The  analysis  is  developed  as  a  continuous-time  Markovian  process  with 
stationary  transition  probabilities.  Probabilities  that  exactly  i  components  have  failed 
by  time  t  are  calculated  based  on  three  levels  of  information:  for  systems  in  unknown 
condition,  for  systems  known  to  be  in  UP-condition,  and  for  systems  whose  condition 
was  not  observed  for  some  amount  of  time.  A  key  part  is  the  investigation  of 
conditional  probabilities  of  i  components  having  failed  by  time  t  for  a  system  known  to 
b".  UP,  the  conditional  limiting  distribution  as  t-»co,  and  relation  to  the  system  failure 
rate.  State  probabilities  for  systems  not  monitored  continuously  for  being  UP  arc 
shown  to  be  bounded  between  those  corresponding  to  systems  that  are  either  observed 
constantly  or  not  at  all. 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 7 

A.  THE  MODEL  OF  THE  3-COMPONENT  STANDBY 

SYSTEM . 7 

B.  IMPORTANT  ASSUMPTIONS  OF  THE  MODEL . 8 

II.  MODELING  THE  3-COMPONENT  STANDBY  SYSTEM  . 10 

A.  TRANSITION  PROBABILITIES . 10 

B.  SYSTEM  SURVIVAL-.  DISTRIBUTION-,  AND 

DENSITY  FUNCTIONS . II 

III.  STATE  PROBABILITIES  IN  THE  3-COMPONENT  SYSTEM  . 12 

A.  STATE  PROBABILITIES  OF  A  SYSTEM  IN  UNKNOWN 

CONDITION . 12 

B.  STATE  PROBABILITIES  FOR  A  SYSTEM  IN  KNOWN 

CONDITION  . 15 

1.  Limiting  State  Probabilities . 16 

2.  Interpretation  of  the  Limiting  Conditional  State 

Probabilities  . 19 

3.  Svstem  Failure  Rate  as  a  Linear  Function  of  the 

Conditional  State  Probabilities  . 20 

C.  STATE  PROBABILITIES  FOUND  IN  A  SYSTEM  NOT 

OBSERVED  CON  TINUOUSLY . 21 

APPENDIX  A:  TRANSITION  PROBABILITIES  FOR  THE  4-STATE 

MODEL . 25 

APPENDIX  B:  SYSTEM  CHARACTERISTICS  . 30 

LIST  OF  REFERENCES . 31 


INITIAL  DISTRIBUTION  LIST 


4 


LIST  OF  TABLES 


I.  LIMITING  STATE  PROBABILITIES  Q.,  X^X^Xg  . 18 

II.  LIMITING  STATE  PROBABILITIES  Q.,  X.,X.  =  X=X  . 19 

i  j  k 

III.  TRANSITION’  PROBABILITIES  STRUCTURE  I  (X^X^X^  . 25 

IV.  TRANSITION  PROBABILITIES  STRUCTURE  2A 

(X0,Xx=X2  =  X)  . 26 

V.  TRANSITION  PROBABILITIES  STRUCTURE  2B 

(X0=X2  =  X,X1)  . 27 

VI.  TRANSITION  PROBABILITIES  STRUCTURE  2C 

(X0  =  X1  =  X,X2)  . 28 

VII.  TRANSITION  PROBABILITIES  STRUCTURE  3 

(X0  =  Xi=X2=X)  . 29 

VIII.  SYSTEM  SURVIVAL-,  DISTRIBUTION-,  DENSITY 

FUNCTIONS  . 30 


5 


a  t'. t-.  j  *_•  »a . 


• .  -  .**  aV 

i_”  \  m.  i  «  *  .  •  k 


L  .  ^  d 


LIST  OF  FIGURES 

Unconditional  State  Probabilities  P^t),  X0:Xx:X2  =  4:2:1  . 13 

Unconditional  State  Probabilities  P  (t),  X0:Xx:X2  =  2:2:2 . 13 

Unconditional  State  Probabilities  P  (t),  X0:Xx:X2  =  1:2:4 . 14 

Conditional  State  Probabilities  Q3(t),  =  4:2:1  . 16 

Conditional  State  Probabilities  Q.(t),  X0:Xx:X  =  2:2:2  . 16 

Conditional  State  Probabilities  Q^t),  XQ:X1:X2  =  1:2:4  . 17 

Composed  State  Probabilities  for  State  1,  XQ=X1  =  X2  =  X . 24 


N"  V  '  ••  .V 


I.  INTRODUCTION 


This  paper  considers  a  non-repairable  system  of  three  components  paying 
attention  to  various  system  structures  that  result  from  the  use  of  non-identical 
components.  The  objective  of  the  analysis  is  to  derive  and  investigate  state 
probabilities1  for  the  system  assuming  different  levels  of  information  available  about 
the  system,  and,  as  far  as  applicable,  with  respect  to  different  arrangements  of  the 
components  within  the  system.  Although  limited  to  a  3-componcnt  system,  the  basic 
results  may  be  extended  to  non-homogeneous  systems  with  more  than  three  elements. 

A  stochastic  model  was  found  to  be  a  convenient  approach  to  the  problem.  It  is 
developed  in  Chapter  II  together  with  the  entire  set  of  transition  probabilities.  That 
chapter  will  also  show  the  system  survival  function,  distribution  function,  and  density 
function  as  functions  that  do  not  depend  on  the  order  of  the  components  within  the 
system.  Chapter  III  will  derive  the  state  probabilities  in  three  steps.  A  subset  of  the 
transition  probabilities  will  immediatedly  determine  the  state  probabilities  of  a  system - 
in  unknown  condition.  Special  emphasis  is  given  to  a  system  which  is  known  to  be  in 
working  condition.  In  this  case  the  resulting  conditional  state  probabilities  will  be 
discussed  in  their  limiting  distribution  and  in  their  relation  to  the  system  failure  rate; 
the  position  of  the  most  reliable  component  in  the  system  will  turn  out  to  be  significant 
for  the  limiting  conditional  state  probabilities.  The  final  step  will  derive  state 
probabilities  for  a  system  which  is  not  monitored  continuously. 

Most  computational  work  is  not  shown  within  the  paper,  however,  the  reader  is 
provided  with  extensive  tables  which  list  the  basic  results  for  all  possible  system 
structures. 

A.  THE  MODEL  OF  THE  3-COMPONENT  STANDBY  SYSTEM 

flic  system  under  consideration  throughout  this  paper  consists  of  three 
components  that  may  be  thought  of  as  an  original  and  two  spare  components.  Only 
one  of  the  components  is  active  at  a  given  point  in  time  and  is  exposed  to  failure. 
When  a  component  fails  it  is  replaced  by  one  of  the  spares,  if  available,  flic 


The  svstem  will  be  said  to  be  in  state  i  if  exactly  i  of  the  components  arc 
DOWN;  thus  the  state  probabilities  correspond  to  the  probability  that  (»,  1.  2,  or  all 
components  have  failed  at  some  time  t  (see  Chapter  II). 
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switch-over  is  assumed  to  cause  no  initial  shock  to  the  new  component  and  to  occur 
immediatedly  after  failure  and  in  a  negligible  amount  of  time. 

The  system  works  (and  is  said  to  be  UP)  as  long  as  at  least  one  of  its 
components  is  operable;  it  fails  (is  DOWN)  if  all  components  have  failed. 

No  repair  facility  is  provided  for  a  broken  component;  thus  it  remains  in  the 
DOWN-condition  throughout  the  mission  once  it  has  failed.  This  implies  that  the 
system  itself  will  fail  with  certainty  in  a  finite  amount  of  time. 

It  will  always  be  assumed  that  the  system  starts  its  useful  life  at  t  =  0  with  all 
components  in  the  UP-condition. 

B.  IMPORTANT  ASSUMPTIONS  OF  THE  MODEL 

The  main  characteristics  of  the  system  arise  from  three  basic  assumptions: 

(a)  Components  in  the  spare  status  cannot  fail.  Therefore  the  failure  rate  of  a 
spare  is  0  (cold  standby)  and  the  lifetime  of  a  component  starts  at  the  moment 
it  is  switched  into  the  active  state. 

(b)  The  active  component  has  an  exponential  life  distribution,  i.e.,  there  is  no 
effect  of  age  on  the  component  failure  rate. 

(c)  All  life  lengths  are  mutually  independent;  i.e.,  the  performance  of  the 
component  currently  active  does  not  depend  on  the  performance  of  its 
predecessor)  s). 

These  assumptions  are  essential  for  the  Markovian  model  used  later  on;  they  establish 
constant  transition  rates  throughout  the  states  of  the  system  with  changes  possible 
only  at  the  time  a  failure  ( =  transition  to  the  next  state)  occurs. 

Note,  however,  that  the  life  lengths  of  components  in  the  active  state  arc  not 
assumed  to  be  distributed  identically.  Thus  the  system  may  consist  of  components  with 
different  expected  lifetimes. 

Whenever  it  is  required  to  distinguish  between  systems  with  respect  to  the 
homogeneity  of  their  components,  then: 

(a)  Structure  1  (\0xk1*Xz)  will  refer  to  a  system  with  three  components  each 
having  a  different  failure  rate  (the  most  general  case); 

(b)  Structure  2  (X..,  X^  =  XR  =  X)  will  denote  a  system  with  two  components  of  the 
same  kind  (having  failure  rate  X)  whereas  the  component  used  in  the  i1*1  place 
has  a  different  failure  rate  X.; 


II.  MODELING  THE  3-COMPONENT  STANDBY  SYSTEM 


The  system  may  be  found  at  any  point  of  time  t  to  be  in  one  out  of  four  possible 
states  {0, 1,2,3}.  It  will  be  said  to  be  in  state  i  if  exactly  i  of  the  components  are 
DOWN;  thus  the  system  will  be  L'P  in  states  {0,1,2}  and  DOWN  in  state  3.  Since  no 
repair  is  provided,  it  is  immediate  that  the  system  after  starting  its  life  at  t  =  0  in  state  0 
will  transit  into  state  1  and  2  and  finally  be  absorbed  in  the  DOWN-state  3. 

The  amount  of  time,  T.,  the  system  spends  in  state  i  before  making  the  transition 
into  state  j  is  a  random  variable,  distributed  exponentially  with  mean  l;Xi,  where  is 
the  failure  rate  of  the  component  active  in  state  i.  By  assumption  1\  is  independent  of 
previous  component  failures  in  the  system. 

The  system,  or  better  yet  the  underlying  process,  may  thus  be  modeled  as  a 
continuous-time  Markovian  process  with  stationary  (homogeneous)  transition 
probabilities  [Refs.  l,2:pp.  26,  234J. 

Moreover,  the  process  constitutes  a  pure  birth  process  with  a  finite  state  space. 
Transitions  are  possible  within  the  four  states  only  from  state  i  to  state  i+  1  and  occur 
at  rate  State  3  is  absorbing  and  allows  no  further  transitions  once  it  has  been 
reached. 

A.  TRANSITION  PROBABILITIES 

By  the  mcmoryless  property  of  the  exponential  sojourn  times  in  a 
time-homogeneous  Markovian  process,  the  transition  probabilities  are  conditioned  only 
on  the  state  of  the  system  at  time  t  =  ()  but  not  on  the  amount  of  time  the  system  has 
already  spent  in  that  state. 

Let  Z(t)  =  0,1, 2, 3  indicate  the  state  of  the  system  at  time  t,  and 
lyt)  =  P[Z(t)  =  j  |  Z(0)  =  i] 

be  the  probability  that  a  process  in  state  i  at  time  t  =  0  will  be  in  state  j  some  time  t 
later. 

Then  the  transition  probabilities  P„(t)  may  be  obtained  by  solving  the 
differential  set  of  the  Kolmogorov  forwards  equations  for  the  pure  birth  process,  as 
done  by  Ross  [Ref.  2:p.  244). 


In  genera!  formulation 

P..(t)  =  e‘V 

ir  ' 

P.  .(t)  =  jVjs  ds 

for  i  =  0, 1 ,2  and  i  +  1  <  j  <  3. 

The  computational  work  is  too  bulky  to  be  presented  within  this  paper,  however, 
the  resulting  transition  probabilities  are  completely  listed  in  Tables  III  .  .  .  VII  (see 
Appendix  A). 

B.  SYSTEM  SURVIVAL-,  DISTRIBUTION-,  AND  DENSITY  FUNCTIONS 

We  will  assume  in  what  follows  that  at  time  t  =  0  the  system  is  new,  that  is  that 
Z(0)  =  0. 

Let  T  be  the  failure  time  of  the  system.  Then  the  probability  that  the  system  will 
survive  a  mission  of  length  t  is  just  the  probability  that  the  system  will  not  transit  into 
the  DOWN-state  before  the  end  of  the  mission 

F(t)  =  P[T>  t]  =  l-P0J(t) 

Furthermore 

F(t)  -  PIT  <  tl  =  P03(t) 
and 

f(t)  =  d/dt  [F(t)|  =  d/dt  [P03(t)| 

The  system  survival  function  F(t),  distribution  function  F(t),  and  the  density 
function  f(t)  arc  summarized  in  Table  VIII  for  the  three  possible  structures  of  the 
system  (see  Appendix  B).  In  cither  case,  the  formulas  obtained  for  the  functions  arc 
completely  symmetric  in  the  indices  of  the  X..  Thus  they  do  not  depend  on  the  order  in 
which  the  various  components  are  used  in  the  system. 


III.  STATE  PROBABILITIES  IN  THE  3-COMPONENT  SYSTEM 


In  the  following,  state  probabilities  for  the  3-component  system  are  derived  for 
three  levels  of  information  about  the  system: 

(a)  State  probabilities  P.(t)  for  a  system  which  are  unconditioned  on  any 
information  about  its  condition  after  the  initial  start-up  at  t  =  0 

(b)  Conditional  state  probabilities  Q^(t)  for  a  system  which  is  known  to  be  in 
UP-condition  at  time  t. 

(c)  Conditional  state  probabilities  P^  (t)  for  a  system  which  was  observed  to  be  in 
UP-condition  at  some  time  s  <  t. 

If  the  system  was  observed  to  be  in  a  specific  state  i  at  some  time  after  start-up 
then  this  information  would  reinitialize  the  Markovian  process  and  either  lead  back  to 
one  of  above  cases  or  to  a  different  model  with  less  components.  Therefore  this  case 
will  not  be  discussed  in  the  following  analysis. 

Most  of  the  expressions  developed  for  the  state  probabilities  in  this  chapter 
contain  several  exponential  terms  (or  a  combination  of  two  or  more  exponential  and 
constant  terms).  Unless  we  restrict  to  simple  cases,  the  attempt  to  calculate  extreme 
points  or  intersection  points  for  the  functions  will  lead  to  equations  which  can  only  be 
solved  using  numerical  procedures.  Therefore,  an  approach  is  used  which  illustrates 
the  functions  in  probability  plots  for  various  parameter  values  and  stresses  special 
features  and  limiting  behaviour  as  t-»so. 

A.  STATE  PROBABILITIES  OF  A  SYSTEM  IN  UNKNOWN  CONDITION 

Within  this  section  it  is  assumed  that  no  information  about  the  state  or  the 
condition  of  the  system  is  available. 

The  probability  that  a  system  will  be  found  in  state  j  at  time  t,  given  it  started  its 
useful  life  at  t  =  0  in  state  0,  is  just  the  probability  of  a  transition  into  state  j  during 
(0,t).  To  simplify  notation,  we  may  define 

Pj(t)  =  P[Z(t)  =  j  |Z(0)  =  0|  =  P0.(t) 

These  state  probabilities  P^(t)  have  been  plotted  in  Figures  3.1  to  3.3  for  three 
failure  rate  combinations  and  over  a  standardized  time  axis  (in  multiples  of  the  system 
MTTF  =  1  1  Xx+  1  X2). 
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Figure  3.1  Unconditional  State  Probabilities  P.(t),  X(J:X1:X2  =  4:2:1. 


Figure  3.2  Unconditional  State  Probabilities  P.(t),  X0:Xx:X2  =  2:2:2. 
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Figure  3.3  Unconditional  State  Probabilities  P^t),  X():X1:X2  =  1:2:4. 

General  information  about  the  state  probabilities  P^(t)  may  be  obtained  by 
considering  their  first  derivatives 

d.'dt  [P.(t)l  =  -X.p.(t)  +  X.^p.^t) 

The  state  probabilities  P.(t)  will  increase  if  the  probability  of  a  transition  into 
state  j  exceeds  the  probability  of  a  transition  out  of  state  j  at  time  t.  Therefore  it  is 
immediate  that  PQ(t)  and  P3(t)  are  monotonously  decreasing,  respectively  increasing 
functions  in  t. 

Pj(t)  and  P,(t)  appear  to  constitute  unimodai  functions"  whose  maxima  can  be 
located  at 

(a)  For  Pj(t): 

t  =  — l—  (lnX1-lnX0)  if  X0<Xx 
and 


"A  proof  for  this  statement  was  not  found  in  the  literature,  however,  all  trials  to 
construct  a  counter  example  have  failed  within  the  work  done  lor  this  study. 
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if  X0  =  Xx  =  X 


(b)  For  P2(t): 

A  closed  form  expression  exists  only  for  a  homogeneous  system: 


if  X0  =  X1  =  X2  =  X 


The  limiting  state  probabilities 

0  if  j  =  0,1,2 

P.  =  lim  P.(t)  = 

J  f  1  if  j  =  3 

confirm  the  expectation  that  the  system  will  finally  be  absorbed  in  the  DOWN-state  3. 

B.  STATE  PROBABILITIES  FOR  A  SYSTEM  IN  KNOWN  CONDITION 

The  following  sections  will  consider  systems  which,  without  revealing  their  exact 
states,  provide  information  enough  to  decide  on  their  overall  condition  —system  UP  or 
DOWN  —  at  time  t. 

If  the  system  is  known  to  be  UP  at  time  t  then  it  may  be  in  either  one  of  the 
states  {0,1,2}  with  probability 


Q3(t)  ■  P[Z(t)  =  j  |  system  UP.  Z(0)  =  0} 


(eqn  3.1) 


j  =  0,1,2 


which  can  be  readily  obtained  from  the  transition  probabilities  in  Tables  III  .  .  .  VII. 
Recalling  that  P0^(t)  =  P3(t),  the  conditional  state  probabilities  may  be  rewritten 


Q.(t)  =  - P.(t) 

J  1  -  P3(t)  J 


j-0,1,2 


and,  for  a  fixed  t,  can  thus  be  obtained  by  multiplying  the  state  probabilities  P3(t),  we 
would  get  without  the  additional  information,  each  by  the  same  factor  c  =  [l  —  P3(t)|'  . 
This  factor  c  is  greater  than  or  equal  to  1  for  all  t  and,  as  a  monotonously  increasing 
function  of  t,  may  get  as  large  as  infinity  as  t-»^o. 

The  probability  plots  for  the  conditional  state  probabilities  Q3(t),  in  Figures  3.4 
to  3.6,  are  based  on  the  same  ratios  in  the  failure  rates  as  the  plots  of  Figures  3.1  to 
3.3  and  may  provide  a  rough  idea  about  the  shape  of  the  functions  Q.(t). 


:X2  =  4:2:1. 


MTTF 


Dbabilities  Q^(t),  XQ:X1:X2  -  2:2:2. 


6 


Figure  3.6  Conditional  State  Probabilities  Q^t),  XQ:X1:X2  =  1:2:4. 

1.  Limiting  State  Probabilities 

The  plot  of  Figure  3.6  suggests  the  existence  of  limiting  conditional  state 
probabilities  that  are  not  necessarily  equal  to  0  or  1,  unlike  the  case  of  the 
unconditional  state  probabilities  P^(t).  State  2  does  not  seem  to  take  over  the  part  of 
an  absorbing  state  in  all  cases,  as  we  might  have  expected. 

To  derive  the  limiting  conditional  state  probabilities 

Q.  =  limQ.(t)  =  Iim  P[Z(t)  =  j  I  system  L'P,  Z(0)  =  OJ 

3  t— »co  3  t— 

for  j  =  0,1,2  it  is  necessary  to  pay  close  attention  to  the  order  in  which  components 
with  different  failure  rates  are  used  in  the  system.  We  also  need  to  specify  which 
component  has  the  minimum  failure  rate.  It  is  therefore  convenient  to  derive  the 
limiting  probabilities  separately  for  the  different  orderings  of  component  failure  rates 
possible  within  the  system. 

a.  Structure  l  (X0*k1*k2) 

In  this  structure  the  failure  rates  of  the  components  are  all  different;  we 
can,  therefore,  specify  a  unique  minimum  failure  rate  It  should  be  recalled  that  the 


indices  of  the  Xi  denote  the  position  of  the  corresponding  component  within  the 
system  (initial  component  or  i1*1  spare);  thus  e.g.,  Xx  =  min(X0,X  ,X  )  will  express  that 
the  first  spare  is  assumed  to  be  the  component  with  the  smallest  failure  rate  used  in  the 
system. 

Using  Equation  3.1  and  Table  III  we  get  the  limiting  probability  matrix 
shown  in  Table  I. 


TABLE  I 

LIMITING  STATE  PROBABILITIES  Q  X0*Xx*X2 

Qo 

— 

Qi 

q2 

X0  =  min(X0,X1,X2) 

X1  =  min(X0,X1,X2) 

X2  =  min(X(J,X1.X2) 

(K-KXK-Ki 

XxX2 

0 

0 

. 

X0(X2-XQ) 

XiX2 

X2-Xx 

K 

o 

K 

K 

i 

_ i 

b.  Structure  2  (X^  X  ^  =  XR  =  X) 

The  limiting  state  probabilities  Q.  may  be  achieved  either  by  calculating 
the  limits  of  the  CU(t),  as  t-»oo,  or  they  may  be  denved  from  the  general  structure 
(Table  I)  using  the  following  rules: 

(a)  If  X.<X,  i.e.,  if  there  is  a  unique  component  associated  with  the  minimum 
failure  rate,  then  use  that  row  of  the  matrix  which  corresponds  to  the  position 
of  the  singular  component  within  the  system 

(b)  If  X<X.,  i.e.,  if  the  minimum  failure  rates  applies  to  two  components,  then 
use  that  row  of  the  matrix  which  corresponds  to  the  position  of  that 
component  with  failure  rate  X  which  is  used  last  within  the  system  (i.e.,  if  two 
rows  in  the  matrix  apply  then  use  the  lower  one). 

The  limiting  state  probabilities  are  listed  in  Table  II. 
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TABLE  II 


LIMITING  STATE  PROBABILITIES  Q.,  X.,X.  =  X  = 

^3’  i’  j  k 

X 

Svstem  Structure 

Qo 

Qi 

Q2 

Case  X.-X-X:  X  <  X 

<v«s 

X0(X-X0) 

^0 

0  0 

X2 

X2 

X 

x0>x 

0 

0 

1 

Case  X-X.-X;  X  <  X 

0 

l  l 

X 

X 

Xx>X 

0 

0 

1 

Case  X-X-X2:  X2<X 

0 

0 

1 

X  >X 

0 

X 

X 

2 

^2 

^2 

c.  Structure  3 

In  this  structure  the  system  consists  of  three  identical  components.  The 
limiting  state  probabilities  may  be  calculated  to  be 

Q0  =  0 
Qj  =  o 
Q2  -  l 

or  they  may  be  looked  up  in  Table  I  by  taking  the  last  row  of  the  matrix  (in 
accordance  to  the  rule  that  points  to  the  last  used  component,  if  more  than  one  has  the 
same  minimum  failure  rate). 

2.  Interpretation  of  the  Limiting  Conditional  State  Probabilities 

One  might  have  expected  to  find  systems  working  on  their  very  last  spare,  if 
they  are  still  alive  after  use  for  a  sufficiently  large  amount  of  time.  However,  the 
results  show,  that  under  certain  arrangements  of  the  component  failure  rates,  we  may 
very  well  find  systems  working  even  on  their  initial  component. 


An  explanation  for  this  system  behaviour  is  most  readily  obtained  by  looking 
at  the  upper  triangular  probability  matrix  in  Table  I  for  a  system  composed  of 
components  all  with  different  failure  rates.  Obviously,  the  component  with  the 
minimum  failure  rate  takes  on  an  important  position  within  the  system.  Given  that  the 
system  is  UP,  the  limiting  conditional  state  distribution  will  put  all  its  probability  mass 
onto  states  corresponding  to  the  most  reliable  component  (main  diagonal  in  the 
matrix)  and  its  possible  replacements.  Thus  functioning  systems  drop  through  initial 
less  reliable  components  at  least  onto  their  most  reliable  ones.  One  will  always  find 
systems  which  remain  working  on  their  most  reliable  component.  Others  will  continue 
to  transit  into  further  states  (unless  the  most  reliable  component  is  also  the  last  one 
available).  If  the  initial  component  is  also  the  most  reliable  one  then  either  one  of  the 
spares  may  be  found  to  carry  the  system  in  the  limiting  distribution;  there  is  no  second 
drop  through'  to  a  spare  with  a  smaller  failure  rate.  The  probabilities  that  a  system 
will  eventually  be  found  in  states  corresponding  to  the  most  reliable  component  and  its 
possible  replacements  depend  only  on  the  failure  rates  associated  with  these 
components. 

In  case  of  the  more  restrictive  structures  2  and  3  we  are  not  always  able  to  . 
define  a  unique  component  to  be  the  most  reliable  one.  Previous  observations, 
however,  remain  valid  if  we  assign  the  attribute  'most  reliable'  to  that  component 
among  identical  ones,  which  is  used  last  within  the  system.  Therefore,  if  more  than  one 
component  shares  the  same  minimum  failure  rate,  we  will  observe  systems  dropping 
through  identical  components  and,  in  the  limiting  probability,  a  system  with  three 
identical  components  will  be  found  working  on  its  last  spare  with  certainty. 

3.  System  Failure  Rate  as  a  Linear  Function  of  the  Conditional  State  Probabilities 

The  system  failure  rate  r(t)  may  be  considered  to  be  the  probability  of  an 
instantaneous  system  failure  at  any  time  t,  given  the  system  survived  up  to  time  t. 

Then 


P[svstem  fails  in  (t,t  + h)|  svstem  I  P  at  tj 
r(t)  =  lim - : : - 

h— *o  h 


=  lim 

h— *0 


& 


P[Z(t  +  h)=3|Z(t)  =  i)  P[Z(t)=i  |  system  UP) 


h 


=  lim 

h-»0 


Pi?(h)Q.(t) 
fsb  h 


=  v  Q  (t)  lim  - 
f?()  1  h— *0 


tyh> 

h 
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Since  the  transition  probabilities 
P.  .(h) 

p. .  —  lim  — u- 
13  h-»o  h 

are  constant  for  a  Markovian  process  with  stationary  transition  probabilities,  the 
system  failure  rate  can  be  expressed  as  a  linear  combination  of  the  conditional  state 
probabilities  Q3(t). 

For  the  3-component  system  under  consideration,  we  have  p()3  =  p13  =  0  and 
p23  =  X2  •  The  system  failure  rate  thus  reduces  to 

r(t)  =  X2Q2(t) 

and  calculates  to: 

(a)  Structure  1 

XlX2(Xl-X2)C'A°t_X0X2^0-X2)e'/’11  + 

(b)  Structure  2 

.  .  2  [(X.-X)t-lle-Xt  +  A1 
rt  1  lXi(X.-2X)  +  XX.(Xi-X)t]c-At  +  k2cKil 


(c)  Structure  3 


r(t)  = 


1'2  X-t 


-h2 


1  +  Xt  +  'iX‘t‘ 


These  system  failure  rates  are  continuous  functions  and,  as  shown  in  general 
by  Barlow  and  Proschan  [Ref.  3:p.  100],  strictly  increasing  over  t.  Thus  the 
3-component  system  wears  out  as  it  ages. 

The  system  failure  rate  r(t)  does  not  depend  on  the  order  in  which  the 
components  are  arranged  within  the  system  (however,  Q.(t)  does).  In  the  limit,  as 
t-»  x  ,  r(t)  always  approaches  the  minimum  failure  rate  found  in  the  system. 

C.  STATE  PROBABILITIES  FOUND  IN  A  SYSTEM  NOT  OBSERVED 
CONTINUOUSLY 

In  the  previous  sections  state  probabilities  were  derived  for  two  special  cases: 

(a)  State  probabilities  P^( t)  for  a  system  that  did  not  reveal  any  information 
about  its  condition  at  time  t. 


(b)  Conditional  state  probabilities  Q^(t)  for  a  system  that  was  known  to  be  in 
L'P-condition  at  time  t. 

We  now  want  to  combine  both  cases  and  calculate  state  probabilities  P^  (t)  for  a 
system  that  was  observed  in  L:P-condition  at  some  time  s  in  the  past  but  whose  present 
condition  (at  time  t)  is  unknowm.  Since  the  system  may  have  failed  during  (s,t)  the 
DOWN-state  3  must  be  included  in  the  analysis. 

As  before,  we  implicitly  assume  that  the  system  was  new  when  started  up  at  t  =  0, 
i.e.,  that  Z(0)  =  0.  By  conditioning 

Pj  (t)  =  P[Z(t)  =  j  |  system  UP  at  s) 

=  y  P[Z(t)  =  j  |  Z(s)=  i]  P|Z(s)=  i  |  sytem  UP  at  s] 

i£j 

=  y  p^t-s)  Qj(s) 


and  for  j  =  0, 1 ,2  by  Equation  3. 1 

p> 

Note:  This  equation  is  also  valid  for  P3  (t),  if  the  sum  is  taken  for  i<j  only  (since 
P[Z(s)=  3  I  system  UP  at  s|  =  0). 

After  calculating  the  sum,  above  equation  reduces  to 


■V"' -  7^  W"  i-»^ 


or,  in  terms  of  state  probabilities  P  (t). 


P.  (t)  =  - P.(t) 

J  l  -  iMs>  J 


j  =  0,1,2 


Again,  the  set  of  conditional  state  probabilities  for  the  UP-states  of  the  system 
may  be  obtained  by  multiplying  each  of  the  state  probabilities  P.(t)  by  some  factor 
c  =[1  ~P,(s)J  .  However,  c  is  a  function  of  s  only  and,  once  s  is  fixed,  it  remains 
constant  over  t  ^  s. 

I  he  time  s.  at  which  the  svstem  was  observed  to  be  UP,  can  be  located  anywhere 


in  the  interval  (<>,t).  If  s  =  0  then  c  =1,  and  P.  (t)  reduces  to  IU(t),  as  expected  (the 
implicit  assumption  Z(0)  =  <)  does  already  include  the  'system  UP'  information).  If  s=t 
then  c  =  c,  and  we  obtain  the  conditional  state  probabilities  Qdt)  investigated  in  the 
last  section. 


Since  Pj(t)  is  monotonously  increasing  over  t,  c  is  always  found  to  be  1  <  c  ^  c. 
Therefore  the  conditional  state  probabilities  P^  (t)  are  bounded  by  P^.(t)  and  Q^t)  such 
that  P  (t)<P  *(t)<Q  (t),j  =  0,1,2. 

Switching  between  P^(t),  P^  (t),  and  Q.(t)  enables  us  to  calculate  the  state 
probabilities  for  different  techniques  used  in  observing  the  condition  of  the  system.  If 
the  condition  of  a  system  is  monitored  continuously  and  found  to  be  UP  over  a  period 
of  time  then  the  state  probabilities  calculated  at  any  time  within  the  interval  will  follow 
the  Q.(t)  curves.  As  soon  as  observations  are  interrupted  the  state  probabilities  are 
determined  by  P^  (t).  If  no  further  observation  is  taken  from  the  system,  then  the 
increasing  uncertainty  about  the  condition  of  the  system  will  force  the  state 
probabilities  associated  with  the  L’P-states  asymptotically  towards  their  lower 
boundaries  P t).  At  the  same  time,  the  DOWN-state  probability 

•  1 

P,  (t)  =  - (I\(t)-P,(s)| 

3  1  -P3(s)  3  3 

increases  from  0  towards  its  upper  bound  P-,(t). 

This  is  illustrated  in  Figure  3.7  for  state  1  of  the  system.  It  is  assumed  that  the 
system  was  initially  monitored  continuously  (in  L'P-condition).  At  time  Sj-JxMTTF’ 
observation  was  discontinued  until,  in  c.  single  observation  at  s2=2xMTTF,  the  system 
was  found  to  be  still  in  UP-condition.  The  probability  plot  for  state  1  (bold  line)  is 
therefore  composed  of  three  partial  plots  using  Q, (t),  for  0<t<Sj,  and  Pj  (t)  with 
s=  Sj  (for  Sj  S  t  <  s,)  and  s=  s,  (for  t  >  s,). 


0 
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MTTF 
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Figure  3.7  Composed  State  Probabilities  for  State  1,  Xq  =  a1-X2  =  X. 


APPENDIX  A 

TRANSITION  PROBABILITIES  FOR  THE  4-STATE  MODEL 


TABLE  IV 


TRANSITION  PROBABILITIES  STRUCTURE  2A  (X„t  Xx  =  X2  =  X) 
P00(D  -  e'V 

Po.w  =  -A— (c-^-c-V) 

k0-K 

P0,(t)  =  {l(X0-X)t-  lle-^  +  e-V) 

{kQ-k) 

„ _  [X0(X0-2X)  +  XX0(X0-X)t|e-^  +  XVV 

03(0  '  O^X? 

Pn(t)  =  e‘Xt 

P12(t)  =  Xte'Xt 

P13(t)  =  l-(l  +  Xt)e'Xt 

P22(t)  =  e‘Xt 
P23(t)  =  l~e'Xt 


P33(t)  =  1 


-^-  y'  ^  mm 


TABLE  VI 


TRANSITION  PROBABILITIES  STRUCTURE  2C  (XQ  =  X1  =  X,  X2) 

P00(t)  -  e'Xl 
P0|(t)  -  Xtc'Xt 

•y2 

P02(t)  =  -r-r-^KV^-ne-^  +  e-^1} 

(A2-A) 

[X2(X2-2X)  +  U2(X2-X)tle-Xt  +  X2c'^2l 


/ 1\ 


TABLE  VII 


TRANSITION  PROBABILITIES  STRUCTURE  3  (X(J=X1=X2=X) 


Poo(')  - 

l>0|(l)  =  Xte-'-' 

P02(t)  =  '/2X2re'Xt 

P03(t)  =  1  -(1  +Xt  +  '/2X2t2)e'^1 

P„(t)  =  e’Xt 
P12(t)  =  Xte’^ 

P13{t)  =  1-{1+Xt)e-Xt 

P22(t)  = 

I\,(t)  =  1  ~e'Xt 


APPENDIX  B 

SYSTEM  CHARACTERISTICS 


TABLE  VIII 

SYSTEM  SURVIVAL-,  DISTRIBUTION-,  DENSITY  FUNCTIONS 


Structure  1  (X^X^Xg): 

Eft)  —  e 

(VMWWl^r e 


-  V  _  ±ohc\t + ha. Za  c*^2t| 


F(t)  =  i  - 


_ _  r  ^l'^2  c-XQt  _  ^0~^'2  C-X1t  +  ^0~^1  C-X.tj 

(^'0'^'l^^'0"^'2^^'l"^2^  ^0  ^1  ^2 


f(‘>  ■ 

Structure  2  (X.,  X.  =  X.,  =  X): 

1  J  K 


E(t)  = 


iX^-IXi  +  XX^Xj-XHIe’^  +  XVV 
1  (X.-X)2 


|X.(X.-2X)  +  XX.(X.-X)t|e’^t  +  XV^it 

1  ‘  (V,V - 


f(t>  *  n\T  IIO-i-lK-  llA'  +  Ai1) 

Structure  3  (XQ  =  Xx  =  X2  =  X): 

F(t)  =  (1  +  Xt  +  'AXnV^ 

I'(t)  =  1  -  (1+Xt  +  '"iXn2^1 
f(t)  =  'iXVc^1 
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